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The magnetically induced Richtmyer-Meshkov instability in a two-component Bose-Einstein condensate is 
investigated. We construct and study analytical models describing the development of the instability at both 
the linear and nonlinear stages. The models indicate new features of the instability: the influence of quantum 
capillary waves and the separation of droplets, which are qualitatively different from the classical case. We per- 
form numerical simulations of the instability in a trapped Bose-Einstein condensate using the Gross-Pitaevskii 
equation and compare the simulation results to the model predictions. 
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I. INTRODUCTION 

Hydrodynamic phenomena in Bose-Einstein condensates 
(BEC) have recently become of interest, in particular, studies 
of the Kelvin-Helmholtz (KH) and the Rayleigh-Taylor (RT) 
instabilities U,-^ have been undertaken. As quantum gases, 
BECs introduce new types of terms into the hydrodynamic 
equations, which may then lead to nonclassical behavior. For 
example, recent studies of shock waves in BECs and quan- 
tum plasmas demonstrated qualitatively different structures of 
the quantum shocks in comparison with the classic gas- and 
plasma dynamics IH-Ql- For the KH and RT instabilities, it 
has been seen that the quantum effects led to an upper cutoff 
for the unstable wavenumbers In addition, BECs al- 

low for wide tunability of the experimental parameters, and 
thus they are excellent systems for exploring large parameter 
regimes and quantitatively testing theoretical predictions. 

Here, we will consider the Richmyer-Meshkov (RM) insta- 
bility, another fundamental hydrodynamic instability in BEC. 
The classical RM instability in gases is closely related to the 
RT instability. The RT instability develops when a heavy gas 
is supported by a light gas in a gravitational field, while the 
RM instability is produced by a pulsed acceleration instead of 
a constant one |8-10]. The acceleration pulse typically oc- 
curs due to a shock wave hitting an interface between two 
gases of different density. In the classical case, the RM and 
RT instabilities demonstrate a qualitatively similar pattern of 
spikes and bubbles with mushroom-structure, though the RM 
instability develops slower, both at the linear and non-linear 
stages. A quantum counterpart of the RT instability in BEC 
was recently demonstrated in 1 3]. Sasaki etal. 1 3] imitated the 
real gravity of the classical case by the gradient of a magnetic 
field in the system of two coupled BECs with different spins, 
e.g. for two hyperfine states of ^^Rb. Numerical simulations 
of the coupled Gross-Pitaevskii equations demonstrated de- 
velopment of the mushroom structure in agreement with the 
classical case. Stabilization of the short-wavelength pertur- 
bations of the quantum RT instability was observed due to 
surface tension between the two BECs. 

The RM instability in the system of a two-component BEC 
subject to a magnetic pulse looks like a natural extension of 
the RT results. Surprisingly, the RM problem of an inter- 
face with surface tension has not been solved yet even in the 
classical case. Only recently, Ref. |11] considered the prob- 



lem but found that the classical methods of the RT-RM theory 
(e.g. the Layzer model) fail when surface tension influences 
the instability. Though brief and qualitative, the discussion 
of Ref. 1 1 1] pointed out that strong changes are expected in 
the instability development at the nonlinear stage because of 
surface tension. In the following work, we show that dra- 
matic changes happen even at the linear stage of the insta- 
bility, which were overlooked in [l lj. The main reason for 
these changes is that the system possesses intrinsic dynami- 
cal properties determined by the capillary waves, and this is 
qualitatively different from the classical RM instability at an 
inert interface. In that sense, the present case resembles the 
RM instability at a flame front studied in Refs. fl2','l3'l, for 
which the intrinsic flame evolution dominates asymptotically 
over the temporary effects of a passing shock. 

In the present paper we investigate the magnetically in- 
duced RM instability in a two-component BEC. We construct 
and study analytical models describing the instability devel- 
opment both at the linear and nonlinear stages. The models 
indicate new features of the instability: influence of quantum 
capillary waves and separation of droplets, which are qualita- 
tively different from the classical case. We perform numeri- 
cal simulations of the instability in a trapped BEC using the 
Gross-Pitaevskii equation and compare the simulation results 
to the model predictions. 



II. ANALYTICAL MODELS 

The dispersion relation for interface waves in a BEC con- 
sisting of two components with different spins in a magnetic 
field gradient B' may be presented as JH] 



fiBB'k/2m 



(1) 



where oj is the perturbation frequency, k = 27r/A is the wave 
number, A is the wavelength of the perturbation, /^s is the 
Bohr magneton, and m is the atomic mass. In the case of 
zero magnetic field, the dispersion relation describes capillary 
waves with frequency 



(2) 



where a is the surface tension and n is the condensate con- 
centration. A positive magnetic gradient in the dispersion re- 



2 



lation leads to the magnetic counterpart of the RT instabil- 
ity with Im(a;) > 0, Re(aj) = for sufficiently long per- 
turbation wavelength. In the present work we are interested 
in the RM instability, which develops in the case of an im- 
pulsive magnetic field jisB' /2m — US{t), where the am- 
plitude U has the dimension of velocity. For this case, the 
dispersion relation Eq. ([T]i should be re-written as a time- 
dependent differential equation with respect to the interface 
position z = f{t) exp{ikx), as 



= [U6{t)k^ul]f. 



(3) 



The original analysis by Richtmyer and Meshkov con- 
sidered the case of zero surface tension cr = 0, which gives 
the capillary wave frequency Uc ~ 0. In the classical geom- 
etry of the RM instability, taking the initial condition f = fo 
att = 0, one can integrate Eq. Q to find linear growth of the 
perturbation, / = Ukfot + fo. Surface tension modifies the 
solution to Eq. (O drastically. Solving Eq. (O for i < we 
find capillary waves 



f = fr expiikx — iuct) + /; exp(iA;a; + iujct), 



(4) 



where labels r and I denote waves propagating to the right 
and to the left, respectively. The capillary waves determine 
the initial conditions for /_ and [df /dt)^ at the interface just 
before the magnetic pulse att — —0. Conditions just after the 
pulse, at i = +0 , are obtained by integrating Eq. Q as 



[/] 



0, [df/dt]t ^ Ukf^, 



(5) 



where [...] designates changes of the respective value during 
the pulse. Still, the solution to Eq. (O after the pulse at t > 
is also a superposition of the capillary waves of Eq. ^ with 
new amplitudes Fr, and Fi. The new amplitudes are related 
to the initial amplitudes using Eq. (|5]l, and are found to be 



F, 



r,l 



Jr,l ± (Jr 



fl)- 



(6) 



Thus, instead of the linear perturbation growth inherent to the 
classical case we obtain energy pumping and redistribution 
in the intrinsic capillary waves of the BEC. In the limit of 
a very strong pulse, Uk/uc ^ 1, we observe the tendency 
of a standing wave to form out of an initial arbitrary wave 
distribution. If we take initial perturbations in the form of a 
standing wave as well with fr = /;* = (l/2)/s exp(ia), fs 
and a being the amplitude and the phase, then Eq. ^ yields 
the new standing wave amplitude Fg as 



Fs=fs 



Uk 



■ cos a 



(7) 



The transformation of magnetic energy to the hydrodynamic 
one is the most effective for the magnetic pulse acting on BEC 
when the standing capillary wave is at it's maximum, a = , 
and this gives 





F, = f,^i + u^kyoji 



(8) 
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Figure 1 : (color online) Detachment of droplets in the RM instability 



In the opposite case, with a — tt/2, the pulse does not 
produce any effect on the interface at all. The combina- 
tion Uk/ujc plays the role of the Weber number. We = 

U^P/u}^ = 2nmU^/ak. 

The linear analysis holds as long as the perturbation am- 
plitudes are sufficiently small, with kFg ^1. At the non- 
linear stage, the classical RM instability leads to bubbles and 
spikes, which resemble qualitatively the nonlinear stage of the 
RT instability. The bubble velocity tends asymptotically to a 
constant for the RT geometry, and it slows down for the RM 
configuration according to 1/kt \ 10]. For the case of a BEC 
we expect a qualitatively different behavior of the instability 
due to two physical mechanisms: 1) The perturbation growth 
is finally stopped and pushed back by surface tension; 2) An 
elongated perturbation finger becomes energetically unstable 
with respect to detachment of droplets. The first mechanism 
may be understood already from the above linear analysis. To 
investigate the second mechanism, we study the system illus- 
trated in Fig. [T] We compare the circumference of the pertur- 
bation finger of amplitude Fi to a similar finger of smaller am- 
plitude F2 plus a circular droplet of diameter d for a fixed total 
surface area of the red part of the figure. Taking, for example, 
a sinusoidal shape of the perturbations, we obtain that droplet 
detachment provides a gain in the surface energy for the initial 
amplitude exceeding the critical value Fi/ D ^ 2.70, which 
corresponds to F2/D « 0.28 and d/D « 0.70. Taking an ax- 
isymmetric counterpart of Fig. 1 we obtain similar estimates 
for the critical ampHtudes Fi/D « 1.42, F2/D « 0.50 with 
the radius of a spherical droplet d/D k, 0.33. Of course, 
such a model provides mainly qualitative understanding of 
the nonlinear stage of the RM instability in BEC, and other 
processes are expected to make the interface dynamics more 
complicated. For example, the secondary KH instability de- 
velops at the interface between two components gliding along 
each other at different velocities Because of the sec- 

ondary KH instability the perturbation fingers may acquire 
a mushroom-shape, and thus spread beyond the limits of the 
channel sketched in Fig. [T] 
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III. NUMERICAL SIMULATIONS 

We take our system to be described by the two coupled 
Gross-Pitaevskii equations (GPEs) 



ih 
ih 



dt 
dt 



2mi 
2m2 



+ V2 (r) + 522 |i/'2 1^ + 512 1 V'l n V': 



(9) 



where 4'j is the macroscopic wavefunction, Vj is the harmonic 
trapping potential on component j, and the interaction param- 
eter is given as 



3D _ ^irh^ajji 



(10) 



with ajj' the s-wave scattering length and rajy the reduced 
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Figure 2: (color online) Density of the lower component, showing the 
capillary waves on the interface at the point of maximum amplitude. 



For our calculations we employ a two dimensional (2D) 
system with a modified interaction parameter, g^^, so that our 
2D simulation can be related to a 3D tight pancake trap. Ori- 
enting our system so the tight direction is in z, our interaction 
parameter is modified as follows I.15il : 



2'Kh'^ajji 

Ijj'Clz 



(11) 



where is the harmonic oscillator length in the z direction. 
We choose our two components to be the two hyperfine states 
of «^Rb, that is iFjup) = |1,1) and |1,-1). These have 
scattering lengths lll4ll an = 022 = 100. 4as and ai2 — 
101. Sas where as is the Bohr radius. 

We begin with a system of 3.2x10^ atoms of ^^Rb, equally 
split between the two hyperfine states. We use the trap geom- 
etry Ux = iOy = 2tt X lOOHz and = 2tt x 5kHz. We solve 
numerically the coupled GP equations to find the ground state 
of the system. The two components are confined to y > 
and y < respectively. Our numerical investigation into the 
RM instability is arranged as follows. First, we induce capil- 
lary waves on the interface between the condensates and study 
their evolution, this is detailed in Sec. IIII A I The standing 
capillary waves are then used as the initial state for the RM 
instability. We trigger the RM instability at the interface, us- 
ing a pulsed magentic field in Sec. IIII Bl and investigate the 
different regimes of the instability. We then make a compari- 
son between our analytical models and numerical results, and 
explore reasons for any discrepancy in Sec. IIII CI 



A. Capillary waves 



another. This results in the amplification of the perturbation, 
and the eventual formation of the RT instability |3]. When the 
waveform on the surface has reached suitable amplitude, we 
switch off the magnetic gradient and allow the capillary wave 
to oscillate freely. 

Surface tension and hence the frequency of oscillations de- 
pends on the density at the interface as cr cx n'^/^, uic oc n^^^, 
with a maximum at the peak density, and reducing to near zero 
at the edges of the system. This introduces finite size effects 
into our system, and for this reason, any quantitative results 
will focus on the system behavior at peak density. We measure 
the frequency of oscillation at this point to be ujc — 246s~^. 
An example of the capillary wave in the system at peak am- 
plitude of oscillation is shown in Fig. |2l While the waves in 
the center of the system are uniform, we observe irregularities 
on the edge of the cloud. The irregularities develop due to the 
variations of the oscillation frequency across the wavefront, 
which lead to a continuous phase shift of the capillary waves 
at the center and on the edges. 

A further effect of the reduced surface tension and inter- 
component repulsion on the edges of the cloud is the 'pass- 
through' of condensate from one side of the system to the 
other Interpenetration of the condensates has been observed 
already in simulations Our numerical results show the 
same effect, though in a reduced form as compared to that in 
1 3]: our choice of initial perturbations allows the interface in- 
stability to develop faster, and so there is less time for conden- 
sate pass-through to occur. We use the configuration shown in 
Fig. |2]as the initial state for studying the RM instability in the 
following sections. 



To form capillary waves at the interface of the two con- 
densates, we introduce a very small sinusoidal perturbation of 
A =9/im to each component along the interface. A steady 
magnetic gradient of B' = 1.78G/cm is added to the sys- 
tem, directed so that the condensates are pushed toward one 



B. RM instability for different pulse strength 



We apply a magnetic pulse of the form B' oc exp (— i^/r^) 
to the system, where r — 0.02 ms controls the duration of the 



4 



pulse. We note that while the field is not a delta function as 
in our analytical model, the pulse duration is much less than 
the characteristic time scale of the capillary wave oscillation, 
uJcT ^ 1- The magnetic "shock velocity" of the analytical 
model is calculated as 



2m 



(12) 



We administer the shock to the initial state shown in Fig. |2l 
The behavior of the system following the pulse depends on 
the size of the magnetic field gradient B' . In agreement with 
our theoretical predictions, we see a wide range of behavioral 
regimes, comprising of both linear and nonlinear dynamics. 
We can further split the nonlinear behavior into three distinct 
regimes: simple nonlinear, nonlinear with detachment, and 
turbulent nonlinear We detail the defining features of these 
regimes in the remainder of this section. 

In the linear phase typical for rather weak magnetic pulses 
with [/ < 3 X lO^'^ms"^, We < 0.70, the pulse results in am- 
plification of the waveform on the interface without qualita- 
tive modifications of its' shape. After a short phase of growth, 
the waveform then reduces in height, and eventually returns to 
capillary oscillations. Throughout the entirety of the evolution 
of the system the wave at the interface remains sinusoidal. 

To demonstrate the simple nonlinear phase, in Fig. [3] we 
show the system at times 6.4ms and 12.7ms after a pulse of 
strength 0.75 x lO^'^ms^^, which corresponds to the Weber 
number We = 4.4. After the pulse, the sinusoidal waveform 
grows at a rate proportional to the pulse amplitude, in agree- 
ment with the analytical model, Eq. ((Sj. As the perturbation 
humps grow beyond the limits of linear dynamics kf <C 1, 
they become distorted, and then eventually 'pinch' in to form 
bubbles, as can be seen in Fig. [3]A. At the later time in Fig. 
|3]B we can see that some of these bubbles (on the edges of 
the cloud) have detached, while others in the middle have re- 
mained connected to the main body of the condensate. The 
difference in behavior happens due to variations of surface 
tension and local Weber number across the interface. Because 
of the lower surface tension the magnetic pulse causes greater 
disturbance on the cloud edges making the flow pattern sim- 
ilar to the classical RM instability lIs UlOll . At the very edges 
of the cloud in Fig. Oa), where surface tension is almost 
zero, the growing perturbations acquire the shape of a half of 
a "mushroom". The clearly visible vortex spiral at the side of 
the "mushroom" cap is a characteristic feature of the classical 
RM instability for gases of comparable density with zero sur- 
face tension. The vortex develops as a result of the secondary 
KH instability of two gases (fluids, condensates) gliding along 
each other because of the primary RM instability. In our case 
the condensates are of the same density, which creates opti- 
mal conditions for the secondary KH instability at the cloud 
edges. Qualitatively different dynamics of the condensate at 
the trap center and at the edges may illustrate the difference 
between the "classical" and quantum RM instability. 

For a stronger magnetic pulse [/ > 8 x lO^^'ms^^, We > 
5.0, the nonlinear phase is accompanied by detachment of 
droplets all over the interface. We present the system at 
times 5.7ms, 8.0ms and 12.7ms after a pulse of strength 




Figure 3: (color online) Cloud density of one component for times 
6.4 and 12.7ms after shock size U = 0.75 x 10"^ms"^ 



We = 8.6 in Fig. 21 which demon- 
strates some new features in the system evolution. In Fig 
|4ta) we find that instead of simple bubble formation we get 
mushroom-shaped caps reminiscent of RT instability. It is in- 
teresting to note the different shape of these mushrooms, as 
compared to those produced due to the RT instability |3]: for 
the RM instability, the mushroom caps have less spiralled den- 
sity under the caps. We conclude that the surface tension has 
a much larger effect for a pulse acceleration, inhibiting the 
growth of the secondary KH instability. The respective phase 
pattern of the system at i = 5.7ms is shown in Fig. |5] We can 
see that there are indeed vortex pairs forming on the caps in 
Fig. Ufa). These vortices then evolve in a complicated man- 
ner, giving rise to results seen for example in Fig. |4lb),(c). 
Another important feature of the phase pattern Fig. |5]is the 
large-scale acoustic oscillations of the trapped condensate as 
a whole (the "breathing" mode). As we demonstrate below in 
Section ITlI CI the acoustic oscillations influence development 
of the RM instability noticeably, leading to quantitative de- 
viation of the numerical results from the analytical model of 
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an incompressible flow without confinement. The distinctive 
feature of the quantum RM instabihty for sufficiently strong 
magnetic pulses is detachment of the mushroom caps from the 
main cloud shown in Fig. EJb). The detachment occurs when 
the distance from peak to trough is approximately twice the 
effective finger width, which correlates reasonably well with 
the nonlinear analytical model, taking into account limitations 
of the model. In the present case the magnetic pulse is not so 
strong, and the droplets are reattached back to the main cloud 
after a short time interval. The process of reattachment is ac- 
companied by formation of small-sized droplets and pockets 
in the condensate components. 

Finally, the instability effects become quite strong for an in- 
tense magnetic pulse leading to turbulent interface dynamics. 
In Fig. |6]we present an example of the highly turbulent system 
12.7ms after a pulse of strength U — 2.1 x 10^'^ms^^ cor- 
responding to the Weber number We — 34.5. This regime is 
characterized by the formation of multiple vortex pairs, and by 
highly complex evolution of the interface. In Fig. |6] one can 
hardly identify separate humps, fingers or mushrooms clearly 
visible in the previous instability regimes. Instead, the inter- 
face takes the form of a mijcing layer resembling the RT cas- 
cade of extreme strength lldilTll . 



C. Comparison with analytical results 

We now compare our simulation results to the analytical 
model. Figure Q shows the analytical predictions of Eq. ([8]l 
for the increase in the hump amplitude because of the mag- 
netic pulse (presented by the solid line) and the respective 
numerical data. In the simulations, we measure the pertur- 
bation amplitude in the very center of the system just before 
the pulse, and the maximal amplitude after the pulse result- 
ing from the RM instability; their ratio is shown in Fig. |2]by 
filled diamonds. Gradient shading in the background indicates 
different regimes of the instability. 

At low pulse strength, in the linear regime, we see good 
agreement between simulation and theory. As the pulse 
strength is increased, in the simple nonlinear regime, the 
simulation results deviate from the theoretical predictions, 
demonstrating noticeably stronger perturbation growth. At 
very high pulse strength, in the nonlinear regime with detach- 
ment of droplets, we see that the theory and simulations once 
again come to agreement. 

The deviation of the numerical data from the theoretical 
predictions is related to the finite size of the system and to the 
acoustic oscillations of the trapped condensate. As we apply a 
magnetic pulse to the system, the condensate gets compressed 
in the direction of the magnetic force (along the y-axis), and it 
expands again after the pulse. Compression produces density 
peak close to the interface, which propagates away in the form 
of a weak shock. When reflected from the system boundaries, 
the shock excites acoustic modes in the condensate. The large- 
scale semi-circles in the phase pattern of Fig. |5ja) correspond 
to the acoustic modes. We demonstrate below that the acoustic 
oscillations pump extra energy in the development of the RM 
instability, thus leading to additional growth of perturbations 
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Figure 4: (color online) Cloud density of one component for times 
5.7, 8.0 and 12.7ms after shock size U = 1.05 x lO^^ms"^ 



in comparison with predictions of the analytical model. 

To study the excitations, we calculate the mean-square 
width in the x direction, {x'^), as a function of time. Evi- 
dently, the acoustic mode produces time oscillations of the 
cloud width. To elucidate the interaction of the acoustic mode 
and the RM instability, we perform two simulation runs: 1) 
First, we apply magnetic pulse to the system with initially flat 
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3.5 



Figure 5: (color online) Phase of portion of a cloud, 5.7ms after 
shock of size U — 1.05 x 10~''ms~^ (a), showing density of a single 
finger (b), and phase of a single finger (c). 




Figure 6: (color online) Cloud density of one component at time 
12.7ms after shock size [/ = 2.1 x 10"^ms"^ 



interface, thus eliminating the RM instability; 2) Second, we 
study a similar pulse acting on a system with perturbed in- 
terface and with the RM instability. Both cases are studied 
in response to a pulse of U — 0.75 x lO^'^ms^^. The re- 
spective acoustic oscillations are presented in Fig. [8] We 
see that for a flat interface (dotted line), the acoustic oscil- 
lations demonstrate a noticeably larger amplitude, also, the 
beat in the frequency of width oscillation seen around 15ms 
signals the presence of more than one mode. In the case of 
a perturbed interface, the oscillation amplitude is markedly 
lower, which indicates that some of the oscillation energy is 
dispersed through the instability on interface. We can also 
compare time dependence for the RM perturbation growth and 
the acoustic oscillations, see Fig. |9] As we can see, the rapid 
initial growth of the interface perturbations (dotted line) is fol- 
lowed by a leveling off of growth at around t = 3ms. The 
subsequent periods of perturbation growth and slowdown cor- 
relate quite well with acoustic oscillations shown by the solid 
line in Fig. |9l In this respect we can see why the growth 
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Figure 7: (color online) Plot of relative fringe growth vs shock mag- 
nitude, U. Analytic prediction (solid line), and numerical simulation 
for total finger growth (filled diamonds) and growth until first level- 
ling off (open squares), see text. 



of RM perturbations in the numerical modeling is noticeably 
larger than that predicted by theory. It is likely that additional 
perturbation growth is provided by the energy of acoustic os- 
cillation pumping the instability. In contrast, the theoretical 
model assumes that the shock wave travels away from the in- 
terface region and never returns. The effect of extra pumping 
may be compared to the classical RM instability in gases un- 
der confinement [H In the geometry of Refs. ifTsifl^ . 
secondary shocks/acoustic waves reflected from the bounding 
walls produce secondary RM instabilities and additional cor- 
rugation of the unstable interface. Pumping of the RM/RT in- 
stabilities by acoustic waves in confined combustion systems 
have been also observed in |20, 21]. To minimize the influence 
of acoustic pumping on the RM growth in our simulations, we 
measure the perturbation amplitude at the instant of the first 
leveling off, which corresponds to the first maximum in the 
acoustic mode. The respective amplitudes shown on Fig. Qby 
open squares, demonstrate a good agreement with the theoret- 
ical predictions of Eq. dHJ for all regimes of the instability. We 
also point out that the difference between two ways of mea- 
suring the final amplitude is minor in the nonlinear regime 
with the detachment of droplets. It is interesting to note that 
droplets typically detach close to the maximum of the acoustic 
mode, during the second contraction of the quadrupole excita- 
tions. Once detached, the droplets are no longer subject to the 
extra instability pumped by the oscillations of the cloud, and 
so they stop growing. At very large times, once the droplets 
reattach back to the main condensate cloud, we note that they 
can begin to grow again. 
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IV. CONCLUSIONS AND OUTLOOK 



t[ms] 



Figure 8: Width of cloud in x direction (parallel to the interface) as 
a function of time for perturbed (solid line) and unperturbed (dotted 
line) clouds. 



Thus, in the present paper we have demonstrated the pos- 
sibility of the RM instability in a quantum system of a two- 
component BEC consisting of hyperfine states of *^Rb with 
spins pointing in the opposite directions. The instability is 
triggered by a magnetic pulse pushing the BEC components 
towards each other. We develop the analytical models of the 
instability at the Unear and nonhnear stages, and solve cou- 
pled GP equations numerically to study the RM instability. 
Both theory and simulations indicate new features of the in- 
stability that are different from the classical case and related 
to quantum surface tension and capillary waves. We obtain 
different regimes of the instability depending on the strength 
of the magnetic pulse: the linear regime, the simple nonlinear 
regime, the nonlinear regime with detachment of droplets and 
the regime of turbulent distortion and mixing of the unstable 
interface. Taking into account the influence of acoustic oscil- 
lations of the trapped condensate, we find good quantitative 
agreement of our analytical model and numerical simulations. 
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Figure 9: Relative growth of fringe as function of time for pulse 
U = 0.75 X 10"^ms"^ (filled circles) and width of the cloud, {x^} 
(solid Une, arbitrary scale in y) 



We point out that experimental realization of the RM insta- 
bility requires a capillary wave of non-negligible amplitude as 
an initial condition. We have shown in this paper how such 
an initial condition can be created in a BEC by applying a 
weaker magnetic field over a longer time, thus making use of 
the RT instability. In a realistic experiment, the initial state 
will be random rather than containing a sinusoidal interface 
perturbation and thus the RT field will excite a spectrum of 
wavelengths rather than a single one; nevertheless, the qual- 
itative features will be unchanged. A single-wavelength in- 
terface wave could conceivably be excited by applying a set 
of locahzed optical potentials along the interface or possibly 
combining an optical potential with a mask to give it the de- 
sired spatial shape. 
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